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Abstract 

The 5*0(4) invariance of the transfer matrix for the one-dimensional Hubbard 
model is clarified from the QISM (quantum inverse scattering method) point of view. 
We demonstrate the 5*0(4) symmetry by means of the fermionic L-operator and the 
fermionic i?-matrix, which satisfy the graded Yang-Baxter relation. The transfor- 
mation law of the fermionic L-operator under the 5*0(4) rotation is identified with 
a kind of gauge transformation, which determines the corresponding transforma- 
tion of the fermionic creation and annihilation operators under the 5*0(4) rotation. 
The transfer matrix is confirmed to be invariant under the 5*0(4) rotation, which 
ensures the 5*0(4) invariance of the conserved currents including the Hamiltonian. 
Furthermore, we show that the representation of the higher conserved currents in 
terms of the Clifford algebra gives manifestly 5*0(4) invariant forms. 
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§1. Introduction 

The one-dimensional (ID) Hubbard model is one of the most important solvable models in 
condensed matter physics. The ground state energy of the ID Hubbard model 

N ^11 
^ = - H I] (cLcm+1. + cl^uCms) + U (""^T " " -), (1-1) 

m=l s=U m=l 

with the periodic boundary condition, 

Cn+U = cL> CN+ls = Cis {S =Ti), (1-2) 

was given by Lieb and Wu ^ by means of the coordinate Bethe ansatz method. Here 
cl^g and Cms are the fermionic creation and annihilation operators with spin s(=t|) at site 
m(= 1, 2, ■ ■ ■ , A^) satisfying the canonical anticommutation relations 

{cL'C„','} = Smm'Sss', {4„cJ^,^,} = {cms,C^'s'} = 0' (1-3) 

and Urns is the number density operator 

Ums = cLcms (s =Ti)- (1-4) 

The parameter U is the coupling constant describing the Coulomb interaction. 

The Hamiltonian ( |1 . 1| ) enjoys two su{2) symmetries One is the spin-s-u(2) generated 

by 

N N ^ N 

=Y1 ^tCmi, S- =J2 cIniCmT, = - (^^T " ^mi) , (1-5) 

m=l m=l m=l 

and the other is charge-SM(2) (ry-pairing su{2)) generated by 

N N 1 TV 

= E i-ircl^cl^, r = E i-ircmicmu = 2 E (^-T + ^mi - 1) . (1.6) 

m=l m=l m=l 



When we assume the periodic boundary condition (|1.2| ), the number of sites should be even 
for the consistency of the charge- sm(2). In this case, the spin-su(2) and the charge-su(2) are 
connected through the partial particle-hole transformation 

Cm\. ^ (-l)'"4i, U ^ -U. (1.7) 

As is well known, these two su{2) are not independent and should be considered as elements of 
a bigger algebra so(4) 0, 

so{A) = su{2)® su{2). (1.8) 

The so(4) symmetry may be the most fundamental property of the ID Hubbard model that 
characterizes the various physical features of the model |^. For example, it was proved by 
EBler et al. |P HT0|] that the Bethe ansatz states of the ID Hubbard model are incomplete and 
have to be complemented by the so(4) symmetry. EBler and Korepin |jll],|12[ showed that the 
elementary excitations of the half- filled band constitute the multiplets of so(4). 
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Several authors have discussed the generahzation of the Lie algebra symmetry so(4) to the 



group symmetry 5*0(4). Following Affleck et al. [|T3|, we introduce the 2x2 matrices 



= f fr" ) . = f .1"' ) , « = 1, ■ • ■ , f . (1.9) 

For convenience, the definition of \E'm in this paper is chosen to be different from the usual one 
P,|l3l. However, they are essentially equivalent. 



The spin-S'f/ (2) transformation can be realized by the left multiplication of an SU (2) matrix 

while the charge-S'f/(2) transformation corresponds to the right multiplication of another SU{2) 
matrix, 

^mC'charge; C'charge ^ SU(2). 

Since the left and the right matrix multiplications are commutative, the transformation 

C^spin^m ^charge (-^"-^0) 

gives the SU{2) x SU (2) transformation among the fermion operators. To put it more precisely, 
the exact group symmetry is 

^0(4) = [SU{2) X SU{2)] /Z2, 

because the choices Cgpin = —1, C'charge = 1 and Ogpin = 1, C'charge = — 1 induce the same 
transformation. The infinitesimal transformation of ( |1.10| ) gives the Lie algebra symmetry 

an 



The integrability of the ID Hubbard model with periodic boundary condition was estab- 
lished by Shastry ||l^-|T6| and Olmedilla et al. |1^,|18|. Shastry introduced a Jordan- Wigner 
transformation 

= {a^ ■ ■ ■ a^_^)a~, Cmi = {(jI ■ ■ ■ (TN){rl ■ ■ ■ t^_i)t~, (1.11) 
which changes the fermionic Hamiltonian ( p..l| ) into an equivalent coupled spin model 

N N U ^ 

= X/ (^m+l^m + ^m^m+l) + X/ i'^m+l^m + '^m'^m+l) + V X/ '-^mTmi (l-^^) 
m=l m=l m=l 

where a and r are two species of the Pauli matrices commuting each other, and and are 
defined by 

1 1 

cr=^ = -fa'^ ± icr^ ) = ±iTy) 

For this equivalent coupled spin model, Shastry constructed the L-operator and the i?-matrix 
(see Appendix), which satisfy the Yang-Baxter relation 

Rl2{0l, 02) [Lm{0l) ® LU02)] = [LrrM ® ^m(^l)] Rl2{01, ^2). (1.13) 

The Yang-Baxter equation for Shastry's i?-matrix was recently proved in |TU|-pi[]. 
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The coupled spin model ( |1.12| ) is also referred to as the ID Hubbard model, since they 
are related through the Jordan- Wigner transformation However, there are differences 



between the coupled spin model ( |1.12| ) and the fermionic Hamiltonian ( |1 . 1| ) . It is well known 
that the periodic boundary condition for the fermion model does not correspond to the periodic 
boundary condition for the equivalent spin model. Moreover, due to the non-locality of the 
Jordan- Wigner transformation (|1 . 1 1|) , the generators of the so(4) symmetry, (|1.5| ) and ( |1.6|) , 
become the non-local in terms of the spin operators a and r. Thus it is more transparent to 
employ the fermionic formulation of the Yang-Baxter relation developed by Olmedilla et al. 



iTj , when we investigate the SO {A) or other symmetries of the ID Hubbard model from the 
QISM (quantum inverse scattering method) point of view. Recently, Gohmann and Murakami 
^ demonstrated that the transfer matrix constructed from the fermionic L-operators has the 
su{2) © su{2) symmetry. The main purpose of this paper is to generalize their result to the 
finite symmetry, namely the 5*0 (4) symmetry corresponding to (|1.1CI|) . 

The plan of this paper is as follows. In section 2, we give a brief summary of the fermionic 
formulation of the QISM for the ID Hubbard model. Some important properties of the fermionic 
i?-matrix are explained. In section 3, we prove the 5*0(4) invariance of the fermionic transfer 
matrix. It is shown that the SO {A) rotation for the fermion operators is related to a kind of 
gauge transformation of the fermionic L-operator. When the number of sites is even, we can 
establish the 5*0(4) symmetry of the transfer matrix under the periodic boundary condition. 
When the number of sites is odd, we have to impose a twisted periodic boundary condition to 
establish the 5*0(4) symmetry of the transfer matrix. In section 4, we discuss the invariance 
of the transfer matrix under the partial particle- hole transformation. In section 5, we give a 
new representation of some higher conserved currents using the Clifford algebra. The 50(4) 
invariance of the conserved currents becomes obvious in this representation. The last section 
is devoted to discussions. 



§2. Graded Yang-Baxter Relation for the ID Hubbard Model 

As a preparation for later sections, we shall summarize the fermionic formulation of the ID 
Hubbard model ||1^|18],^. The fermionic L-operator is 

fmi{0) 



-i/mT(^)cLi 



— fm^{0)Cml 



-hJ 



le cl^i^cl^]^ 



is CjTi^ Cjnj. 
(0) 

9m'\{(^)Cml 



\ 



-e^9m]{0)9mi{0) J 



(2.1) 



where 



fms{0) 
9ms{0) 



sin 6 — {sin 6 — i cos 6} n„ 
cos 9 — {cos + i sin 9} n„ 



The parameter h is related to the spectral parameter 9 and the Coulomb coupling constant U 
through the relation 

sinh2/;, U 

= — ■ 2.2 

sin 2^ 4 ^ ' 



We express by the Grassmann (graded) direct product 



P(1)=P(4) = 0, P(2)=P(3) = 1. 
The fermionic L-operator satisfies the graded Yang-Baxter relation 



under the constraint of the spectral parameter 

sinh 2hi sinh 2h2 
sin 201 sin 262 



The explicit form of 7^12(6*1, 6*2) is |17 



7^l2(e 

/ a+ 















V 



where 



h) = 

e 



-ib+ 
















e 






-i5^ 














d+ 




i/ 



-i/ 



-c+ 






ib- 
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-i/ 



d- 




— 




i/ 













e 






ib- 









ib- 





e 
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1' 







i/ 



—c" 



d- 




-i/ 

















e 



\b- 








-c+ 



-i/ 



i/ 



rf+ 




















e 





COS 




3^(^i-^2)|l±tanh(/;,i-/i2) 
sin(^i - Q2) cos(^i - ^2) |l ± tanh(/ii - h 

{COS 
1 ±tanh(/ii + h2) 
cos I 

sin2(^i - ^2) |l ±tanh(/ii + h2y^ 
y sin( 

1 ± tanh /ii - /i2 , 

COS^fc*! + O2) 

, , , sinful — 62) 

l±tanh /ii + /i2 ■ )n\nv 
sm(6'i + O2) 

cos(^i - ^2) ^ ^ sin(gi - ^2) 

cosh(/;,i + h2) 















-i6^ 


e 





(2.3) 
(2.4) 

(2.5) 



\ 


















(2.6) 



cosh(/ii — h2 



(2.7) 



4 



The second equalities for the Boltzmann weights 6^ and are vahd due to the constraints 



(|2.5|) . There are useful relations among the Boltzmann weights IjT^ 

rf± = a± + c±, = a+a- + b+b-, f = b+b~ + c+c-, 

a+c+ + a-c- = + (r)^ a+c" + a-c+ = 26+6", 

a+ + a~ + 0-^ + c- = 2. (2.8) 

For convenience, we introduce an equivalent fermionic i?-matrix 

nu{ei,e2) = Vi2ni2{0i,e2), (2.9) 

where Vu is the graded permutation 

V^,,f,s = (2.10) 
In terms of 7^.12(^1,^2) ( P-9D , the graded Yang-Baxter relation ( p.4| ) is cast into 



7^l2(^l,^^2) (/:^(^^i)®/) (/®/:^(^2)) = (/f i:™(^2)) (i:^(^l)®/)7^l2(^l,^^2). 

Here / is the 4x4 identity matrix. 

The fundamental properties of the fermionic i?-matrix 7^12(^1, 62) are summarized as follows 

il- 

(1) Regularity (Initial condition): 

ni2{9o,9o)=Vi2. (2.11) 

(2) Graded Yang-Baxter equation: 

ni2i0i, ^^2)7^l3(^^l, ^3)7^23(^^2, ^s) = 7^23(^^2, ^3)7^l3(^l, e3)nuiei, ^2), (2.12) 

under the constraints 

sinh2/ii sinh2/i2 sinh2/i3 U 
sin 2^1 ~ sin 2^2 ~ sin 2^3 ~ T' 

(3) Unitarity: 

7ei2(^l, W2l(^2,^l) =P(^1,^2) /, (2.13) 

where 

7^2l(^2,^l) =Pl27^12(^2,^l)Pl2, 

and 

p(9i, 62) = cos2(0i - 62) (cos2(0i - 62) - tanh2(/ii - /i2) cos2(^i + 62) 



Since the non-zero elements of the i?-matrix ( p.9|) are even with respect to the parity P{a) 
(E3), i.e., 

P{a) + P{f3) + P{a) + P{I3') = (mod 2) for 7^„^^,,^, (^1, ^2) 7^ 0, 
the graded Yang-Baxter equation (|2.12|) can be expressed in terms of the matrix elements as 



W7''(^2,^3)7^„y^."V(^l,^3)7^„",'^.','(^l,^2)(-)^(^■■^(^(")+^^^^ (2.14) 



Here the summentions are taken over the repeated indices. 



In our previous work [0], we found two important relations of the fermionic i?-matrix with 
constant matrices M and V . The first relation is the symmetry of the fermionic i?-matrix 



0, 



(2.15) 



(2.16) 



7^l2(^l, ^2), M®M 

s 

where the general form of the matrix M is given by 

M 

with the condition 

AM = M11M44 - M41M14 = M22M33 - M23M32. 

We call the matrix M symmetry matrix. For simplicity, we assume AM = 1 throughout the 
paper. 

The second relation is 











Mi4 


\ 





M22 


M23 








M32 


M33 







V M41 








M44 





v®v 



7ei2(^i,^2;f/) 

where the general form of the matrix V is given by 



V(^V 



n 



02-, -U), 



(2.17) 



V 



( 

V21 

V 



V12 



V^2 



^13 




V^43 



\ 

V^34 

/ 



~ ^13V42 — V21 V34 — V31 V2 



'13 1^42 



21 V34 



-^31 y24:- 



(2.18) 



In the relation (|2.17| ), the ^/-dependence of the fermionic i?-matrix is explicitly written. The 
coupling constant of the fermionic i?-matrix in the RHS is —U, or equivalently, hi — *• —hi 
and h2 —>■ — /i2- The matrix V is related to the partial particle- hole transformation (|1.7| ). The 
constant matrices M and V play an important role in the consideration of the symmetry of 
the transfer matrix for the ID Hubbard model (see §3 and §4). We remark that the symmetry 
matrix of Shastry's i?-matrix is not of the form (|2.16|) (see Appendix). This gives one of 
the reasons why the fermionic formulation employed in this paper is more appropriate for the 
investigation of the ID Hubbard model (|1.1|) . 

The monodromy matrix is defined as the ordered product of the fermionic L-operators 



m=l 



T{9) = n i^nm = 1^N{e) ■■■Ci{9). (2.19) 

we have the global relation for the mon- 

(2.20) 



From the (local) graded Yang-Baxter relation 
odromy matrix 

7^l2(^l,^2)[m)®2^(^2)] = \T{e2)®T{ei)]ni2{ei,e2 
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or equivalently 
where 



7^l2(^l,^2)T(^l)f (^2) 



--T (^^2)T(^l)7^l2(^l,^2), 



(2.21) 



T {e) = T{ei)®i, T{e) 



/®T(^2). 

s S 

Define the (fermionic) transfer matrix by 

strKT{e) = tr {((x" ® a') KT{e)} , (2.22) 
where the constant matrix K assumes the form ( |2.16| ) and determines the boundary condition 



24| . In particular, K = I corresponds to the periodic boundary condition (|1.2| ). Then from 
the global graded Yang-Baxter relation ( p.20|) , we can prove that the transfer matrix ( |2.22|) 
constitutes a commuting family 



[strKT{ei),stiKT{e2)] = 0, 



(2.23) 



which proves the integrability of the ID Hubbard model with the (twisted) periodic boundary 
condition. 



§3. SO (4) Symmetry of the Fermionic Transfer Matrix 



We shall discuss the 50(4) symmetry of the fermionic transfer matrix ( |2.22| ). Let us consider 
the following transformation of the fermionic L-operator 

CUd) = M''CUO)M, (3.1) 

where the constant matrices M and M have the form of the symmetry matrix ( |2.16| ), i.e., 



M 



and 



M 



( Mn 








M22 





M32 


V M41 





-M41 




/ Mn 








M22 





M32 


V M41 





t-M4i 


Mi4 



Mi4 

M23 

M33 

M44 J 



'33 



M23M32 



(3.2) 



Mi4 

M23 

M33 

M44 J 



'32 



(3.3) 



Since the matrices M and M are the symmetry matrices, the transformed L-operator Cm{d) 



(|3.1|) also satisfies the graded Yang-Baxter relation with the same fermionic i?-matrix, 

n^2{ei.e2)[trn{e{)®c^{e2)] = [i:™(^2)®>c„(^i)fe(^i,^2). (3.4) 
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We now look for a special transformation of (|3.1| ), which satisfies 

C-mifi'i Cms) — C-mifi'i Cms)- (3-5) 

Here we explicitly write the dependence of the fermionic L-operator on the fermion operators. 
The fermion operators Cms and Cms are assumed to be connected through the transformation 
law (|1.10|) . We discovered that the relation (|3.5| ) is satisfied when the matrices M and M meet 
the following conditions 



M44 


= M*n, 


M41 


= -M*4, M33 


= M2*2, M32 


= -m;. 


Mil 


= Mn, 


M44 


= M44, Mi4 = 


-Mi4, M41 


= -M41 


M22 


= M22, 


M33 


= M33, M23 = 


M23, M32 = 


M32. 



(3.6) 

The condition ( p.2|) now becomes 

|Mn|' + |Mi4r = IM22I' + |M23r = 1- (3.7) 
It is useful to introduce the submatrices of the matrices M and M as 

_( Mu Mu\ ^ _ 

Then the conditions ( |3.6| ) and ( p.7| ) are equivalent to the r( 

Mcharge = O'^M'chargeO'^, Mgpin = Mgpin, Mcharge, Mgpin G S'f/(2) . (3.9) 

The corresponding transformation law of the fermion operators is 



M22 


M23 


M32 


M33 


M22 


M23 


M32 


M33 


ions 





4i i5mT M*2 -M23 \ / 4i iCmT \ ( Mn M14 

i^T j I M22 j 1^ i4^ Cmi ) \ -Ml, Ml, 



(3.10) 

Hereafter, we implicitly assume the the conditions ( p.9| ) for the matrices M and M. Then the 
transformation (|3.1|) is not a gauge transformation in a strict sense, because M ^ M (particu- 
larly Mcharge 7^ Mcharge)- So we try assiguiug the different transformation laws to the L-operators 
for odd sites and even sites as 

£2„-l(^^) = M-^C2n-l{e)M, C2n{e) = M'^C2n{e)M. (3.11) 

The corresponding transformation law of the fermion operators on odd sites is, of course, given 
by the formula ( p.lO| ), 



cL-U iC2n-lT \ = 1 ^2*2 -M23 \ / 4n-U iC2n-lT \ ( M14 
icL-lT ^2n-U / V ^^23 M22 ; V icL-lT ^2n-U / V ~^14 



(3.12) 



Since the matrices M and M are related by the exchange M14 <-h> — M14, the transformation 
law for even sites reads 

( 4ni i52nT \ M*2 -Af23 \ / 4ni iC2nT \ -M14 ^ 

V i4nT £2ni ; ^ MI3 M22 J \ icl^ C2„i ^ ^ ' 
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Multiplying the Pauli matrix cr^ from the right, we obtain 



IC 



2nt 



22 



23 




-lC2nT 



Mil Ml 



14 



-M*4 M*i 



(3.13) 



Recalling the definition of the 2x2 matrices \E'r 
laws (|312D and ( g7[3l) as 



( |1.9D , we can summarize the transformation 



"^m = M^p^^^'^Mcharge, 



(3.14) 



which exactly coincides (|1.1CI| ) with the correspondences Cgpin = M^pin and Ocharge = ^^charge 



As we have explained in §1, the transformation ( p.l4|) is the SO (4) rotation in the space of 
the fermion operators. Therefore we can conclude that a kind of gauge transformation ( |3.11|) 
induces the SO (4) rotations for the fermion operators. We call the transformation ( p.ll|) the 
5*0(4) rotation for the fermionic L-operator ( |2.1| ). Note that the canonical anticommutation 
relation (|1.3|) is preserved under the transformation (|3.14j) 



^ms 1 s J 



mm ss 1 



0. 



Let us consider the 50(4) invariance of the fermionic transfer matrix ( |2.22| ). First we assume 
that is even and impose the periodic boundary condition. The local 50(4) rotation for the 
fermionic L-operators (|3.11 ) induces the 50(4) rotation for the monodromy matrix ( 2.19|) 



f{e) = n ^m{0) = M-'T{e)M. 



(3.15) 



m=l 



Since the relation 

str{X(^)M} = stT{MX{e)} (3.16) 
holds, the transfer matrix ( p. 221) is invariant under the periodic boundary condition [K = I) 

stYf{9; Cms) = strT(6'; c^s), (3.17) 

where we write the fermion operators explicitly. In the relation ( 3.16| ), X{6) is any 4x4 
matrix, which may depend on the fermion operators. On the other hand, the transfer matrix 
strT(^; 

Cms) can be expressed as 

strf(^;c„,) = strT(^^;g^,), (3.18) 
due to the property (p.5|). Combining ( p.lTj) and ( |3.18|) , we establish 

strT(6'; Cms) = strT(6'; Cms)- (3.19) 



The relation ( p.l9|) shows that the fermionic transfer matrix is invariant under the 50(4) 
rotation for the fermion operators ( p.l4| ). It indicates that all the higher conserved currents, 
which are embedded in the transfer matrix, also have the 50(4) symmetry (see §5). 

Now we shall write the transformation (|3.15|) in terms of the submatrices Mcharge and Mgpin 
(|3.8|) . We introduce the following convenient notation for the monodromy matrix 



T{9) 



( Dn{e) Cn{e) Cu{e) Du{e) \ 

Br^{e) An(^) A^2{e) B^^^e) 

B2i{9) A2i{9) A22{9) B22{9) 

V D2,{9) C2i{9) C22{9) D22{9) J 



9 



where we regard A{e) = {Aij{9)), B{0) = {Bij{0)), C{9) = iQj{9)) and D{e) = {Dij{e)) as 
2x2 matrices. Then the transformation ( p.l5| ) can be expressed in terms of 2 x 2 matrices 
A{e),--- ,D{e) as 

A{e) = M-LA(^)M3pi,, 

C{d) = M-pl„C(^^)M,harge, 

b{e) = M,l,Lg,Z}(^)Mcharge. (3.20) 

Because the transformed monodromy matrix also satisfies the graded Yang-Baxter relation with 
the fermionic i?-matrix, 



-^12(^1, 62) 



f(^2)®T(^l 



R 



O2), 



the associative algebra defined by the graded Yang-Baxter relation should be invariant under 
the transformation ( p.20|) . From ( |3.20 ), we notice an interesting fact that the submatrix A{6) is 
transformed by the spin-S't/(2) rotation and D{6) is transformed by the charge-5't/(2) rotation. 
We believe this property plays a significant role in the application of the algebraic Bethe ansatz 
for the ID Hubbard model [P^. 

Next we consider the case of odd. The monodromy matrix transforms as 



f{e) = M-^T{9)M. 



(3.21) 



In this case, we have to twist the periodic boundary condition to make the transfer matrix 
5*0(4) invariant. The condition for the matrix K in the transfer matrix is 



KM = MK. 



For example, 



K 



( i 











\ 





1 
















1 















— i 


/ 



(3.22) 



(3.23) 



solves the condition ( |3.22| ). From the formula in ||2^, we can see that the choice ( |3.23D corre- 
sponds to the twisted boundary condition 



N+ll — Cn+11 



-iCiT, 
-icu- 



(3.24) 



Assuming ( |3.23[ ) and ( ^.24| ), we can prove the 5*0(4) invariance of the transfer matrix for 
odd as 

stiKTiO; Cms) = stiKTiO; Cms), (3.25) 



in a similar way to the even case. 
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§4. Partial Particle-Hole Transformation of the Fermionic Transfer 
Matrix 



In , the transformation law of the fermionic L-operator corresponding to the partial particle- 
hole transformation ( p.. 71 ) was found. We shall discuss the transformation law in connection 
with the relation ( |2.17| ). Consider the following transformations of the fermionic L-operators 



where 



V 



f -1 0\ 

10 

1 

V -10/ 



V 



f 1 \ 

10 

-1 

V -1 / 



(4.1) 



(4.2) 



Since the constant matrices V and V are of the form ( p.l8|) , we have the following graded 
Yang-Baxter relation with the transformed L-operator (^.1|) 



ni2i0i,O2;-u)[cuOuU)®jCUO2;U)] = [jCU02;U)^cuOuU)]ni2iOue2-,-u), (4.3) 



where we write the [/-dependence explicitly. The graded Yang-Baxter relation ([4.3| ) implies 
that the transformed L-operators Cm{0', U) are related to the L-operators with the coupling 
constant —U. In fact the following relations hold, 



C'2n-l{d] C2n-ls, U) 
^2n\d'i C2ns) U) 



'^C,2n-l{d'i C2n~ls, —U), 
^^2n{(^', C2ns) 



(4.4) 



where 



C2n-\\ 
C2n-\ 



C2n-l\-, C2n-li 



-a 



t 

2n-li; 



C2nh 



C2nl 



-'2ni- 



(4.5) 



The transformation ( [4.5| ) is nothing but the partial particle-hole transformation ( |1.7| ). Therefore 
we call ( [4.1|) the partial particle- hole transformation of the fermionic L-operator (|2.1|) . 

It is quite interesting to note that the transformation ( [4.5| ) can be written in terms of the 
2x2 matrix ([T^ ) as 

= -^L, (4.6) 

where f denotes the hermitian conjugation. Moreover, taking the hermitian conjugation of 
(|3A|), we find 

= M-h?,,g,^™M3pi„, (4.7) 

which shows that the spin-S'f/(2) Mspin and the charge-Sf/ (2) Mcharge are exchanged after the 
partial particle-hole transformation. 

We are ready to verify the invariance of the transfer matrix of the ID Hubbard model 
under the partial particle-hole transformation. First we assume that N is even and impose the 
periodic boundary condition. Then the partial particle-hole transformation of the monodromy 
matrix induced by (|4.1|) is 

f{e- c™„ U) = V-^T{e- Cms, U)V. (4.8) 
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From the relations 

strr(^; Cms. U) = str {V''T{e; c^s, U)V} = -strT(^; c^s, U), 

and 

strT(6'; Cms, U) = i^strT(6'; Cms, -U), 



we obtain 

strT(e; Cms, U) = -i^strT(^; Cms, -U). (4.9) 



The last identity ( |4.9| ) proves the invariance of the fermionic transfer matrix (up to sign) under 
the the partial particle-hole transformation ( |1.7|) . Note a relation 



str {X{d)V] = -str {VX{9)} , (4.10) 

which should be compared with ( |3.16|) . 
We have a similar relation for odd, 

str {KT{e; Cms, U)} = -i^-^str {KT{e; Cms, -U)} . (4.11) 

where K is given by (|3.23| ). In the derivation of ( [4.11| ), we have used the relation 

KV-^ = iV-^K. 

The formula ( ^.11| ) means that the fermionic transfer matrix for odd is also invariant under 
the partial particle-hole transformation (|1.7| ) when we assume the twisted periodic boundary 
condition ( p.24| ). 



§5. SO (4) Symmetry of the Higher Conserved Currents 

In §3, we have shown the 5*0(4) symmetry of the transfer matrix, which means that all the 
conserved currents of the ID Hubbard model also have the 5*0(4) symmetry. As will be seen, 
the 50(4) symmetry of the conserved currents can be manifestly read out in terms of the 
Clifford algebra. Hereafter, for simplicity of explanation, we assume the number of sites is 
always even and impose the periodic boundary condition. 

Although the graded Yang-Baxter relation ensures the existence of infinitely many higher 
conserved currents in involution, it is not an easy task to obtain their explicit forms from the 
transfer matrix (see |1^). To construct the higher conserved currents, we often use the boost 



operator p7|,|28|, which recursively produces the higher conserved currents. However in the 



case of the ID Hubbard model, the boost operator does not exist [29| and we have to resort to 
a more direct computation. 



The first higher conserved current of the ID Hubbard model was found by Shastry ||14|,|1| 

as 

AT 

— i {^(^m+2s^ms (^ms^rn+2s^ 

m=l s=ti 

N 

—iU ^ ^ (cm+igCms — CmgCm+ls^ (^m+l,-s + nm,-s — 1) • (5.1) 
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Subsequently, some higher conserved currents were obtained in a similar fashion p9|-pT|, 

N 

/(3) 



5Z 5Z ('^m+3s'^ms ~^ ^\ns^rn+ 
m=l s=U 

N r 

I (^Cjy^_^_igCfn-ls "I" Cm-1< 
m=ls=tl I 



gC-m+ls J I ""m+l,— s ~l~ "'m,—s ~l~ "■m— 1,— s 



n 



m+ls 



-\\n 



'"m,—s 



N ( 
m=l I 



27 ' 



r(4) 



m=l s=ti 

^ / IW 1 

N 

m=l s=tJ. 

Af r m+3 ^ 

m=ls=TI I 



27 v''"^''"' 2 



k=m 
m+2 



A:=m— 1 
m+2 



fc=m— 1 



N ( 

m=ls=ti I 



'TIS 



n. 



m+l,—s 



+ in. 



N 



™=i s=TJ. 



(5.2) 



(5.3) 



These currents are embedded in the fermionic transfer matrix ( |2.22| ) and should be 5*0(4) 
invariant from the result in the previous section. One can confirm the SO {4) invariance of 



13 



these currents using the transformation law of the fermion operators ( p.l4| ). In the following, 
we present a different approach; we shall rewrite these currents in manifestly 5*0(4) invariant 
forms by use of the Clifford algebra. 



Define (j 



and 



where n = 1, 



^ 2n-l — 



■2n-lt 



,4) by 

+ C2n-lt? r 



2n-l 



(4 



2n-l 



■2n-U + C2n-U, T 



2n-l 



2n-lT C2„-l| I , 



i (4n-U - C2„_U 



(5.4) 



rL = i( 



2n 



i (C2n 



2n 



2n 



C2nt + '^2nT; 
4ni + C2ni, 



(5.5) 



, Y- Then the operators satisfy the defining relations of the Clifford algebra 



{^^^^t} = 25,,5'^^ j,k 



(5.6) 

In terms of F^, the SO (4) rotation for the fermion operators (|3.14| ) can be expressed simply as 



j 



^G'^^'Fj, G = {G"^) e SO{A). 



b=l 



The relation between the matrices G and M in §4 is explicitly given by 

f eo ei ^3 -6 \ 

^0 —^2 —^3 
^^3 ^2 ^0 
V 6 ^3 -^1 ^0 / 





Co 


Ci 


-Ca 


-C2 


\ 




-Ci 


Co 


C2 


-C3 






C3 


-C2 


Co 


-Ci 




V 


C2 


C3 


Ci 


Co 





where C,i and rji are real numbers given by 

Co = Re(Mn), 6 = Im(Mn) 
Co = Re(M22), Ci = Im(M22) 



6 

C2 
3 



Re(Mi4), 
Re(M23), 



^3 = Im(Mi4), 
Ca = Im(M23), 



EeJ = Ec? = i. 

j=0 j=0 



(5.7) 



(5.8) 



Clearly the matrix G^^^ corresponds to the charge-S't/(2) and the matrix G'^'^^ corresponds to 
the spin-S'f/ (2). It is an interesting exercise to confirm that the matrices G*-^-* and G*^^^ commute 
each other 

Define the operator F^ by 



■p5 -pi ■p2 T-\3 T-\4 \ ^ ^ -f^a -pfe -pc -pd 

h - h h h h - 1^ ^'^bcd ij i-j ij i-j. 



(5.9) 



a,--,d=l 
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The operator has the following properties 
It is clear that the operators such as 



(j + k), 



Epa pa pE 



(5.10) 



(5.11) 



a=l 



are invariant under the 5*0(4) rotation (|5.7| ). We make use of this fact to rewrite the conserved 
currents. The Hamiltonian 7i = 1*^^^ in terms of the operators V°- and is ^] 



(5.12) 



where 



u 



Here and hereafter, we do not mind the difference of an overall factor. The formula (|5.12|) gives 
a manifestly 5*0(4) invariant representation of the Hamiltonian. 

In the same way, we express the higher conserved currents in terms of the operators and 



T] as 



/(2) 
j(3) 



E E ri+2ri + « E E(-i)'r •+ir," (r,^i - r 
EE(-i)^r^\3r- 

--E|r5 + Er-+2r^(r^2-r5) 

I \ ^ -pa -pa -pfe pbp5 ^"P^ "P^ pap6p5\ 

i+2J- j+iJ- j+i^ j+i - 2 ^+1 i+i i J iy 



(5.13) 



j{4) 



+«^EE(-iyr^+irir5+iFf 
EEr-+4r," 

+«E(-i)4 E {r-+3r- {v%, - r?) - f^^.f^ {t%, - Ff)} 

j L a 

I /"P'^ P*^ P^ p6p5 _|_ pa pa p6 p6p5 

+ \J- i+sJ- i+2J^ jJ- j+2 + J- j+sJ- i+lJ^ i+lJ^ jJ- j+1 

pa p6 pa p6p5 _i_ pa p6 pap6p5l I 

J- i+2J- j+2J- i+lJ- j+2 ^ ^ i+lJ- jJ- jJ- i/ r 

I 2 \ ^ J \ ^ pa pap5 p5 \ ^ pa pa p6 p6p5 / p5 p 

2^^ j 2^ ^ j+l^ j+l^ j+1 V i+2 ^ 

+^'EE(-i)^r-+iri(r,^i-r 



(5.14) 



(5.15) 
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Since the terms that constitute ( |5.13| )- (|5.15|) are of the form ( |5.11| ), we can see that the higher 
conserved currents j(2)^/(3) g^^d I^^^ are also manifestly 5*0(4) invariant. Note that the con- 
straints a 7^ 6 in the summations do not break the 5*0(4) symmetry. For example, we can 
write 



E"pa -pa -pb "p6"p5 \ ^ -pa -pa -pb -pb-pb \ ^ -pa -pa-pb 

J- j+2'- j+lJ- i+lJ- j'- j+1 -2^^ 3-^2^ i+lJ^ i+lJ^ 3^ j+1 J- j+2J- 

ai^b a,b a 

Both terms in the RHS of ( |5.16|) are clearly 5*0(4) invariant. 

The infinitesimal generators of the 5*0(4) rotations (|5.7| ) are given by p3| 



(5.16) 



Q"" = -Q 



ba 



In fact the generator ( 5.17| ) fulfills the defining relation of the Lie algebra so(4) 

'Q^\ Qcdj _ _• (^^bcQad _ ^acQbd _ ^bdQac ^ ^adgbc^ 



We also have a relation 



(5.17) 



(5.18) 



(5.19) 



which is nothing but the infinitesimal transformation of (|5.7|) . Using (|5.19|) , one can confirm 
the commutativity 



Qab^ j(n) 



0, 



n 



(5.20) 



which shows the Lie algebra so(4) symmetry of the conserved currents Actually the 

generators of the spin-sM(2) ( |1.5| ) and the charge-s-u(2) ( |1.6| ) are related to Q"''' as 
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5^ = (Ql2 _ g34 



(5.21) 
(5.22) 



where we introduced 5^, 3^,7]^^ and rj"^ through the relations 

= S"" ± iS^ 7]^ = 7]"" ± ir]^. 
For the Clifford algebra ( ^.61) , the partial particle- hole transformation (|1.7| ) corresponds to 



(5.23) 



Note that the transformation ( |5.23| ) exchanges the spin-su(2) ( |5.21| ) and the charge-su(2) ( |5.22D . 

4 

As for the conserved currents, the transformation ( |5.23| ) preserves the operators like ^ F" F^, 



a=l 



but changes the sign of F^. From the explicit formulas ( |5.12|) - (|5.15|) , one can immediately 
find that the conserved currnets I^'^\n 
transformation 



1, ■ ■ ■ , 4) are invariant under the partial particle-hole 



u 



-u. 



(5.24) 



This is consistent with the result in 54. 
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§6. Discussions 



We have investigated the 5*0(4) symmetry of the ID Hubbard model from the QISM point 
of view. Our approach is based on the fermionic formulation of the Yang-Baxter relation 
for the ID Hubbard model found by Olmedilla et al. It consists of the fermionic -R- 



matrix and the fermionic L-operator. We have discovered the transformation law ( p.llj ) of 
the fermionic L-operator under the SO {4) rotation. It is a kind of gauge transformation and 
induces the transformation of the monodromy matrix. The result is a fundamental property 
of the associative algebra defined by the fermionic i?-matrix of the ID Hubbard model. We 
believe that the property will also play an important role in the algebraic Bethe ansatz for 
the ID Hubbard model, which was recently explored by Ramos and Martins [EBl. We like 
to emphasize the advantage of the fermionic formulation of the Yang-Baxter relation. It is 
difficult to discuss the 5*0(4) symmetry of the Hubbard model through Shastry's i?- matrix and 
the related transfer matrix. 

The 5*0(4) invariance of the transfer matrix ensures the 5*0(4) invariance of the conserved 
currents. We have demonstrated the 5*0(4) symmetry of the higher conserved currents employ- 
ing the Clifford algebra, which corresponds to the spinor representation of the rotation group. 
It should be interesting to explore a representation of the fermionic L-operator itself in terms 
of the Clifford algebra. 

On the infinite lattice, the Lie algebra so(4) = 5^(2) © su{2) symmetry of the ID Hubbard 
model is extended to the Yangian Y(so(4)) = Y{su{2)) (BY{su{2)) symmetry 



Y(so(4)),/« 







as was discovered by Uglov and Korepin |^. The generators of Y(so(4)) can be expressed in 
terms of the Clifford algebra as follows 



Q 



(0) 
ab 



QlS = -iE(-i)'(r°+ir; + r°r 



b 

j+1 



\j>k k>j / c^a,t 



By use of the fundamental properties of the Clifford algebra (|5.6|) , (|5^) and (|5.10|) , we have 
confirmed that the higher conserved currents I^'^\n = 2,3,4) also have the Yangian Y(so(4)) 
symmetry, i.e.. 



nio) Tin) 



nW An) 
^ab 5 



0, 



n 



A. 



All the conserved currents of the ID Hubbard model on the infinite lattice are conjectured 
to have the Y(so(4)) symmetry. In fact Murakami and Gohmann [^] recently showed the 
existence of an infinite number of the conserved currents which have the Yangian symmetry 
on the infinite lattice. However, one of the two Y(s'u(2)) that constitute Y(so(4)) drops out 
p5| . It seems to be difficult to prove the full Y(so(4)) symmetry of the conserved currents 
simultaneously in their method. 
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Appendix : Symmetry Matrix of Shastry's i?-Matrix 

The L-operator for the coupled spin model (|1.12| ) |[l^-[17| is expressed as 



/ 
















^ ^ m' m 








^ ^ m m 


e-'^p-{e)q^{e) 




V 


^ ^ m' m 









(A.l) 



where 



pUo) 



- (cos 6 + sin 6) ± - (cos 6 — sin 6) cr^, 

- (cos 6* + sin 6*) ± - (cos 9 — sin 9) r^. 



(A.2) 



Here the coupling constant h is related to the spectral parameter 9 through the formula ( |2.2| ). 
The i?-matrix -Ri2(6'i, ^2), which satisfies the Yang-Baxter relation ( p..l3| ) with the L-operator 
( |A.1| ), is connected to the fermionic i?-matrix (|2.6| ) through the formula |T^ 



Rl2{9l,92) = Wr2'ni2{9l,92)Wi2, 

where W12 is a diagonal 16 x 16 matrix 



diag(l,l,-i, -i, -i, -i, 1,1, 



-!)• 



(A.3) 



(A.4) 



We consider the symmetry matrix M of Shastry's i?-matrix -Ri2(6'i, 6*2) which is defined to 
be a constant matrix satisfying 



R^2i9^,92),M®M 



0. 



(A.5) 



Here the matrix elements of M are assumed to be commuting numbers. One might suppose the 
symmetry matrix of Shastry's i?-matrix is identical to that of the fermionic i?-matrix (|2.16|) . 
However, surprisingly enough, we notice they take different forms. In fact, solving the defining 
equation for the symmetry matrix ( |A.5| ) directly, we find that the foUowings are the symmetry 
matrix of Shastry's i?-matrix. 



M 



f a \ 
0/300 

7 

V 5 y 



/aOOO\ /00 0a\ 



0/30 
7 

V 5 y 



0/300 
7 

V 5 y 



/ a \ 

0/30 

7 

V 5 y 



(A.6) 
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where a, (3, 7 and 5 are c-numbers obeying 



a5 = Pj. (A.7) 

We ignore the difference of overall factors of the matrices. Then each matrix ([A .61) depends 
only on two parameters. The result means that Shastry's i?-matrix does not reflect the SO (4) 
symmetry of the fermionic Hamiltonian Therefore the SO (4) symmetry of the transfer 

matrix that we explored in this paper may not be discussed if we use Shastry's i?-matrix and 
L-operator (|A.1|). 
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